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FOREWORD 


This  report  describes  a  portion  of  the  results  obtained  on  NASA  Grant 
NSG  3044.  This  work  was  done  under  subcontract  to  the  University  of 
Illinois,  Urbana,  with  Prof.  S.S.  Wang  as  the  Principal  Investigator.  The 
prime  grantee  was  the  Massachusetts  Institute  of  Technology,  with  Prof. 

F.J.  McGarry  as  the  Principal  Investigator  and  Dr.  J.F.  Mandell  as  a 
major  participant.  The  NASA  -  LeRC  Project  Manager  was  Dr.  C.C.  Chamis . 

Efforts  in  this  project  are  primarily  directed  towards  the  development 
of  finite  element  analyses  for  the  study  of  flaw  growth  and  fracture  of 
fiber  composites.  This  report  presents  a  method  of  analysis  for  adhesive 
or  interlaminar  cracks  which  propagate  in  the  interface,  rather  than  co¬ 
hesively  in  the  adhesive  or  interlayer.  The  latter  case  was  treated  in  an 
interim  report  NASA  CR-135248. 
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ABSTRACT 


A  study  on  the  elastic  behavior  of  interface  cracks  in  adhesively 
bonded  lap  shear  joints  is  presented.  The  problem  is  investigated  by 
using  a  recently  developed  method  of  analysis  based  on  conservation  laws 
in  elasticity  for  nonhomo geneous  solids  and  fundamental  relationships  in 
fracture  mechanics  of  dissimilar  materials.  The  formulation  leads  to  a 
pair  of  linear  algebraic  equations  in  mixed-mode  stress  intensity  factors. 
Singular  crack-tip  stress  intensity  solutions  are  determined  directly  by 
information  extracted  from  the  far  field.  Stress  intensity  factors  and 
associated  energy  release  rates  are  obtained  for  various  cases  of  interest. 
Fundamental  nature  of  the  interfacial  flaw  behavior  in  lap-shear  adhesive 
joints  is  examined  in  detail. 


1.  INTRODUCTION 


While  adhesively  bonded  joints  have  been  employed  in  many  secondary 
structural  components,  the  effective  use  of  adhesive  joining  technology 
in  primary  load-bearing  structural  members  is  still  in  its  infancy.  Ap¬ 
plications  of  adhesive  joints  will  not  reach  their  full  potential  until 
certain  critical  problems  are  solved.  The  major  factors  responsible  for 
this  situation  are  the  complex  failure  modes  and  mechanics  presented  in 
the  joints.  Because  of  many  geometric,  material  and  fabrication  variables 
involved,  the  fracture  problem  of  adhesive  joints  is  extremely  complicated. 
Current  knowledge  of  the  joint  failure  behavior  is  yet  very  limited.  Ana¬ 
lyses  and  design  criteria  for  flaw  resistance  of  adhesively  bonded  joints 
are  seriously  deficient. 

Interfacial  cracking,  also  called  debonding,  is  one  of  the  most  fre¬ 
quently  encountered  modes  of  failure.  It  frequently  occurs  at  geometric 
boundaries  such  as  edges  and  re-entry  corners  due  to  inherent  stress  con¬ 
centrations,  or  results  from  faulty  joining  in  fabrication  such  as  in¬ 
complete  wetting  between  adherend  and  adhesive.  Debonding  is  also  found 
in  adhesive  joints  subjected  to  combined  high  temperature  and  moisture 
absorption.  Progressive  reduction  of  stiffness,  exposure  of  the  interior 
of  the  composite  to  environmental  attack,  and  final  disintegration  of  the 
structure  are  frequently  observed  to  occur  by  the  presence  and  growth  of 
interfacial  cracks.  The  lap-shear  joint  shown  in  Fig.  1  is  considered  in 
this  study,  because  it  is  one  of  the  most  widely  used  joint  configuration 
in  structural  applications.  Understanding  the  fundamental  nature  of  inter¬ 
facial  cracks  is  of  utmost  importance  in  the  reliability  and  safety  design 
of  the  adhesive  joints  for  primary  load-bearing  structures.  In  this  paper 
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analytical  study  on  the  interfacial  crack  behavior  in  the  adhesively  bonded 
lap-shear  joints  is  presented. 

The  interface  crack  problem  was  apparently  first  studied  by  Williams 
[1].  Elasticity  theories  for  cracks  of  this  nature  were  attempted  by  many 
investigators,  and  solutions  characteristically  involve  oscillatory  stress 
singularity  [2-7].  A  major  concern  for  the  oscillation  of  near-field 
stresses  is  that  they  lead  to  physical  absurdity  of  crack  surface  inter¬ 
penetration  or  overlapping,  as  pointed  out  by  England  [2]  and  Malyshev,  et 
al.  [8].  This  implies  that  the  solutions  for  the  interface  crack  problem 
are  physically  inadmissible.  The  unsatisfactory  aspects  of  the  oscillatory 
stress  singularity  were  discussed  by  several  researchers,  for  example. 

Refs.  [2,8-10].  An  alternate  model  of  a  closed  crack  tip  with  friction¬ 
less  contact  was  first  proposed  by  Comninou  [9]  to  correct  these  effects. 
The  stress  singularity  of  the  model  is  different  from  that  of  the  oscilla¬ 
tory  one  with  a  finite  normal  tensile  traction  ahead  of  the  crack  tip.  The 
crack  tip  in  a  nominal  tensile  field  has  an  extremely  small  contact  area 
in  comparison  to  the  size  of  the  crack.  Recently,  Achenbach,  et  al.  [10] 
and  Keer,  et  al.  [11]  proposed  another  models  which  include  the  crack  face 
closure,  relative  slip  conditions  at  the  interface,  and  different  defini¬ 
tions  of  stress  intensity  factors.  In  view  of  the  complexities  aforemen¬ 
tioned,  the  interfacial  crack  problem  in  adhesively  bonded  joints  is  ob¬ 
viously  rather  difficult.  The  complex  structural  geometry,  the  presence 
of  the  multiphase  material  system,  and  the  unknown  multiaxial  stress  state 
acting  on  the  crack  introduce  additional  complications.  The  very  local 
nature  of  the  interface  crack  singularity  as  shown  in  [7,12]  is  noted  to 
introduce  the  thickness  of  the  thin  adhesive  layer  as  a  characteristic 
dimension  in  the  problem.  Any  simplifications,  which  fail  to  consider  the 
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crucial  role  of  the  thin  adhesive  layer,  would  lead  to  a  severe  drawback 
in  studying  the  full  nature  of  the  adhesive  joint  fracture.  The  importance 
and  complexities  of  interface  debonding  in  adhesive  joints  have  long  been 
recognized  [13-15] .  However,  research  progress  on  this  kind  of  real-life 
defect  has  been  relatively  slow. 

In  this  paper,  a  study  on  the  interfacial  crack  behavior  in  adhesively 
bonded  lap-shear  joints  is  presented.  By  using  the  recently  developed 
conservation  laws  of  elasticity  for  nonhomogeneous  solids  [16-20]  and 
fundamental  relationships  in  fracture  mechanics  of  interface  cracks,  the 
problem  is  reduced  to  a  pair  of  linear  algebraic  equations,  and  stress 
intensity  solutions  can  be  determined  directly  by  information  extracted 
from  the  far  field.  This  feature  makes  the  current  approach  particularly 
suitable  and  attractive.  Solutions  are  obtained  for  adhesively  bonded 
lap-shear  joints  with  various  material  systems  and  geometric  configurations 
of  the  crack  and  the  joint.  The  method  of  analysis  is  of  practical  use  in 
the  design  and  analysis  of  adhesive  joint  fracture.  The  fundamental  nature 
and  unique  features  of  the  interface  flaw  behavior  in  adhesively  bonded 
joints  are  revealed. 
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2 .  FORMULATION 


The  Eshelby-Rice  conservation  law  for  a  homogeneous  solid  in  a  plane 
elasticity  problem  [16-19]  has  recently  been  extended  to  a  solid  composed 
of  two  dissimilar  materials  [20]  with  the  following  form: 


J.{S}  « 

l 


(W  n. 


-  a 


jk 


uj,i)ds ' 


( [W  ]  6 


i2 


])ds  =  0, 


(1) 


where  W  is  the  strain  energy  density;  cr  ,  the  stress  tensor;  u.,  the  dis- 

J  k  J 


placement  vector;  n^,  the  outward  unit  normal  of  an  arbitrary  closed  contour 
S  which  encloses  a  portion  of  the  continuum,  and  l  is  a  portion  of  the  inter¬ 
face  bounded  by  S  as  shown  in  Fig.  2.  The  bracket  [  ]  in  Eq.  1  denotes  the 
jump  of  a  function  across  the  interface  =  0,  i.e., 


[W]  =  W(xr0+)  -  W(xr0“), 


(2) 


[u.  .]  =  u  (x  ,0+)  -  u  (x  0-). 
J  j1-  J  1  J  y  -1- 


(3) 


Continuity  conditions  of  displacements  and  interlaminar  stresses  across 
x^  =  0  require 

and  [o  2^  =  0-  (4a, b) 


[u. ]  =  ° 


The  component  of  the  conservation  integral  in  Eq .  1  can  be  simplified  to 


J1(S}  = 


(W  nl  -  ajk  “k  Uj,l)dS  =  °* 


(5) 


which  is  the  same  as  the  analogous  result  for  a  single  phase  material. 

For  a  solid  containing  an  interface  crack  between  two  dissimilar  materials 
as  shown  in  Fig.  3,  in  Eq.  5  along  a  path  T  has  the  standard  J-integral  form 
in  homogeneous  fracture  mechanics,  i.e,, 
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j  =  j  {r}  = 


(6) 


3u. 

(Wdx2  -Tii^ds)’ 
r 

where  F  is  an  arbitrary  path  surrounding  the  crack  tip,  provided  that  the 
crack  surfaces  are  free  from  traction  and  the  interface  is  a  straight  line. 

[Path  independent  integrals  (similar  to  Rice's  J-integral  of  course  with  appro¬ 
priate  modifications)  for  elastic  media  with  spatially  varying  moduli  were  noted 
by  Atkinson  [21]  through  an  energy -momentum  tensor  formalism.]  The  J  is  shown 
to  relate  to  the  energy  release  rate G of  an  interface  crack  in  a  usual  manner,  i.e. 

J  =  G  =  G].  +  G  ,  (7) 

where  Gj  and  G^  are  the  energy  release  rates  associated  with  the  mode  I  and 
mode  II  stress  intensity  factors.  The  interface  crack-tip  stress  intensity 
factors,  KI  and  K^,  are  defined  in  a  complex  form  as 

K  -  i  K  =  2/27  e87r  lim  zi+lg  ^(z),  (8) 

z->o 

where  $  (z)  is  a  complex  potential  in  the  well-known  Kolosov-Muskhelishvili 
formulation  [3,22], and  8  is  a  bimaterial  constant  given  in  Ref.  [3]. 

It  is  noted  that  K  and  defined  in  Eq .  8  for  an  interfacial  crack 
in  dissimilar  media  are  different  from  those  for  a  crack  in  a  homogeneous 
solid.  Thus,  the  K  and  may  not  possess  the  usual  significance  and 
physical  interpretation  as  in  the  cohesive  (or  homogeneous)  fracture.  While 
the  overall  stress  intensity  factor  [3,23]  may  be  used  to  express  the  maximum 
amplitude  of  stresses  and  to  correlate  crack  extension,  accurate  description 
of  the  oscillatory  singular  stress  field  in  the  neighborhood  of  the  crack  tip 
requires  detailed  knowledge  of  individual  stress  intensity  factors. 

Because  of  the  complexity  of  the  problem,  numerical  procedures  are  necessary 
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for  obtaining  accurate  solutions  for  interface  crack  problems.  For  example, 
Lin  and  Mar  [24]  used  a  hybrid-stress  finite  element  method  to  achieve  this. 
Hong  and  Stern  [23]  employed  a  contour  integral  method  based  on  Betti's 
Reciprocal  work  theorem  and,  recently,  Smelser  [25]  used  crack-flank  dis¬ 
placement  data  provided  by  numerical  solutions  to  yield  the  interface  crack- 
tip  stress  intensities.  Each  of  these  numerical  approaches  requires  its  own 
computational  scheme  to  handle  the  problem  and  has  given  satisfactory  results. 
In  this  section,  an  alternate  and  convenient  method  of  analysis  is  proposed  for 
determining  the  interfacial  K^.  and  in  bonded  joints. 

Using  the  stress  field  established  by  Rice  et  al.  [3],  one  can  readily 
show  that  the  J-integral  in  Eq.  7  is  related  to  the  stress  intensity  factors 
for  an  interfacial  crack  by 

J  *  j,  nir1®! +  Kix)  <9> 

i=l  i 


for  a  plane  strain  case.  The  coefficients  in  Eq .  9  are  replaced  by  l/[4p. 

(1  +  v^) ]  for  a  plane  stress  condition.  The  J-integral  alone  does  not  provide 
adequate  information  for  determining  individual  values  of  K^.  and  for  an 
inherently,  mixed-mode  interfacial  crack.  Further  development  of  Eq.  9  based 
on  the  introduction  of  known  auxiliary  solutions  for  the  crack  problem  can 
remove  this  difficulty. 

Consider  two  independent  equilibrium  states  with  associated  field  variables 
denoted  by  superscripts  1  and  2  for  the  elastically  deformed  bimaterial 
body.  Superposition  of  the  two  equilibrium  states  leads  to  another  equilibrium 
state,  0.  The  J-integral  for  the  superimposed  state  can  be  shown  to  have  the 
form 


J  (0)  =  j*-1)  +  j('2-)  +  m^1,2\ 


(10) 


7 


where  another  conservation  integral  with  the  form 


M 


(1,2) 


r 

U>2)  dx,  - 

l 

a  (2)  ,  (1) 
(1)  i  (2)  x 

■  \ 

ds 

4 

l  3x^  '  J‘i  3x^ 

r* 

i 

(ID 


The  W^1,2^  in  Eq.  11  is  the  mutual  potential  energy  density  of  the  bimaterial 


body,  defined  by 


W(1>2)  =  Cijki  4)]  4]l  =  Cijk£  4 )]  4)1 


(12) 


Recalling  the  J-K  relationship  in  Eq.  9  for  the  superimposed  state  and  com¬ 
paring  it  with  Eq.  10,  one  can  obtain 


M 


(1,2)  =  2a 


(1)  (2)  (1)  (2) 

KI  KI  +  KII  KII 


(13) 


where 


2  i 

l  ~~a  -  v.). 
i-i  4ui 


(14) 


The  M-integrals  in  Eqs .  11  and  13  deals  with  interaction  terms  only, and  are 
used  directly  in  solving  the  interface  crack  problem.  The  M-integral  is 
clearly  related  to  the  details  of  the  stresses  and  deformation  at  the  crack 
tip  (i.e.,  K^.  and  in  Eq.  13),  but  yet  may  be  evaluated  in  the  far  field 
(i.e.,  the  integral  in  Eq.  11),  where  such  a  calculation  can  be  carried  out 
with  greater  accuracy  and  convenience  than  near  the  crack  tip.  It  is  also 
noted  that,  while  K^.  and  characterize  the  controversial  near-field  oscil¬ 
latory  singular  stresses,  the  energy  release  rate  G,  and,  perhaps,  and  G^, 
are  quantities  well  defined  and  can  be  evaluated  conveniently  mathematically 
and  physically. 

'The  M-integral  used  here  and  elsewhere  in  this  paper  is  defined  by  Eq .  11, 
and  should  not  be  confused  with  those  in  Refs.  [17]  and  [18]. 
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3 .  SOLUTION  PROCEDURE 


Equation  11  together  with  Eq.  13  provides  sufficient  information  for 
determining  the  stress  intensity  solutions  for  a  mixed-mode  interface  crack 
problem,  when  known  auxiliary  solutions  are  introduced.  Denote  the  first 
auxiliary  solution  by  a  superscript  2a  for  a  crack  between  two  dissimilar 
materials  subjected  to  mode-I  loading  only  with 


and 


Equation  13  can  be  simplified  to 


K 


(2a) 

II 


0  . 


(15) 


M 


(1,2a) 


2a  K 


(1) 

I  » 


(16) 


where  M^1,2a^  has  the  same  form  as  Eq.  11  with  the  subscript  2  being  replaced 
by  2a.  T  f  ^  and  uf^  in  Eq .  11  can  be  determined  along  a  properly  selected 
integration  path  V  in  the  far  field  by  any  convenient  method  such  as  the 
commonly  used  finite  element  analysis.  For  a  plane  crack  problem  with  the 
loading  of  Eq.  15,  T^2^  and  u^2a^  are  derivable  by  the  Kolosov-Muskhelishvili 
formulation,  and  the  auxiliary  solution  has  the  following  form: 


u(2a) 


9;  p 


l’ 


(17) 


with 


T  (2a) 


a  (2a) 

ij 


n . 
J 


> 


a(2a) 

ij 


(£n  r. 


9, 


yl’  vl*  p2’  +  °^1')  ’ 


(18) 

(19) 


Exact  forms  of  the  functions  f^^  and  may  be  found  in  [3,4]. 

The  second  auxiliary  solution,  denoted  by  the  superscript  2b,  for  the 
bimaterial  solid  under  a  pure-mode  II  loading  such  that 
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and 


(20) 


K 


(2b)  = 
II 


can  be  found  by  the  same  formulation.  The  conservation  integral  is,  there¬ 
fore,  expressed  by 


Ma>2b> 


2a  K 


(1) 

II  * 


(21) 


where  jt^’2b^  has  a  form  similar  to  that  of  Eq.  11. 

(2a)  (2a) 

It  is  important  to  note  that  the  auxiliary  solutions,  u)  ,  T )  and 
uf2b\  T_£2b\  are  independent  of  the  particular  boundary-value  problem 
being  posed.  Therefore,  they  may  be  determined  independently  by  any  con¬ 
venient  analytical  method  once  for  all.  In  solving  for  and  ,  the 

integrals,  M^1,2^  and  M^1,2b\  must  be  evaluated  accurately  and  explicitly. 
For  a  given  crack  geometry,  loading  condition  and  bimaterial  constant,  this 
can  be  achieved  by  integrating  Eqs .  16  and  21  along  a  properly  selected 
contour  in  the  far  field  so  as  to  avoid  crack-tip  complications.  In 
conjunction  with  the  auxiliary  solutions  determined,  a  numerical  method, 
using  a  conventional  finite-element  approach,  is  currently  employed  to 

calculate  T^  and  u.  .  The  M-integrals  are  then  formally  evaluated  by 
ii 

using  the  second-order  Gaussian  quadrature  along  a  contour  T  passing  through 
Gaussian  stations  of  each  element  (Fig.  4). 
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4.  NUMERICAL  RESULTS  AND  DISCUSSION 


The  solution  procedure  aforementioned  is  programmed  for  studying  general 
two-dimensional  interface  crack  problems.  Evaluation  of  the  conservation 
integrals  is  conducted  in  conjunction  with  a  conventional  finite  element 
method  by  using  eight-node  isoparametric  elements.  Solutions  are  obtained 
for  lap-shear  joints  with  various  adherends ,  adhesives,  and  geometric  con¬ 
figurations.  Accuracy  and  convergence  of  the  results  are  affected  by  several 
unusual  features  of  the  problem  and  the  method  of  analysis  due  to  the  singular 
nature  of  the  crack  and  inherent  approximation  involved  in  the  numerical 
evaluation  of  the  conservation  integrals.  Assessments  of  solution  accuracy 
are  made  by  examining  relevant  problems  for  which  unquestionably  correct 
and  exact  solutions  are  available  in  the  literature.  Excellent  agreement 
is  obtained  between  current  results  and  reference  solutions.  Details  of  the 
study  are  reported  elsewhere  [26],  Current  results  have  an  accuracy  within 
approximately  three  percent  deviation  from  reference  solutions  based  on  the 
optimum  mesh  and  the  integration  contour  presently  selected.  The  primary 
objectives  of  this  section  are  to  determine  stress  intensity  solutions, 
which  characterize  local  stresses  and  deformation,  and  to  examine  effects 
of  material  and  geometric  variables  on  the  crack-tip  response.  Of  particular 
interest  are  the  complex  failure  modes  and  energy  release  rates  associated 
with  the  interface  crack. 

4.1  Interfacial  crack  and  Joint  Configuration 

The  lap-shear  joint  considered  in  Fig.  1  is  composed  of  two  high¬ 
stiffness  and  high-strength  adherends  bonded  by  a  thin  adhesive  layer.  The 
upper  and  lower  adherends  and  the  adhesive  layer  are  assumed  to  have  uniform 
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thicknesses  of  t^,  t^  and  ,  respectively.  The  overlap  region  has  a 
dimension  L .  A  crack  of  length  a  is  located  along  the  interface  of  the 
upper  adherend  and  the  adhesive  near  the  traction-free  edge.  Except  for 
the  interfacial  crack,  perfect  bonding  is  assumed  everywhere.  The  two 
adherends  are  made  of  materials  with  elastic  constants  E^,  and  E^, 
v and  the  thin  adhesive  layer  has  properties  of  E^  and  •  The  adhe¬ 
sive  and  adherends  are  assumed  to  be  linear,  elastic.  Studies  of  related 
cases  such  as  the  center-of-bond  cohesive  crack  and  the  eccentric  crack 
problems  were  reported  elsewhere  [27-28]. 

4.2  Effects  of  Adhesive  Properties  on  Interfacial  Crack-Tip  Response 
Geometric  variables  in  the  problem  are  given  as  the  following: 

6  =  45°,  L  =  0.5  in., 

t^  =  t^  =  0.05  in.,  t^  =  0.005  in.,  (22) 

a  =  2.5  t2  =  0.0125  in. 

The  upper  and  lower  adherends  are  made  of  the  same  aluminum  with  elastic 
constants,  E^  =  E^  =  10  x  106  psi,  and  =  0.33  .  Effects  of 

different  kinds  of  adhesive  on  the  crack-tip  response  are  examined  by 
considering  various  values  of  .  For  the  purpose  of  generality,  adhes¬ 

ives  with  a  broad  range  of  modulus  values  covering  three  decades  on  a 
logarithmic  scale  are  studied.  The  cases  with  =  20  v  40  are  typical 

for  aluminum/epoxy  systems.  Higher  values  correspond  to  the  joints 

with  less  rigid  adhesives  or  subjected  to  a  "hot  and  wet"  environment. 

Crack-tip  stress  intensity  solutions  and  associated  energy  release  rates 
are  obtained  as  functions  of  the  modulus  ratio  shown  in. Figs.  5  and  6. 

Failure  modes  in  the  joints  are  observed  clearly.  The  interfacial  crack 
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experiences  a  mixed-mode  fracture  even  under  uniaxial  loading.  Both 
and  decrease  rapidly  as  the  adhesive  modulus  decreases,  and  the 

difference  between  and  increases  with  E^/E^.  The  mode-I  stress 

intensity  factor  is  higher  than  that  of  the  mode-II  in  the  entire  range 
studied,  except  for  the  cases  of  E^/E ^  ~  1.  This  suggests  that  the  domi¬ 
nant  mode  of  failure  is  opening-mode.  In  the  case  of  E^/E  =  1  ,  the  crack 
is  cohesive,  and  the  nature  of  the  crack-tip  singularity  changes.  The  values 
of  Kj  and  in  this  case  are  found  to  be  consistent  with  the  results 

in  [27],  The  total  energy  release  rate  and  the  energy  release  rate  of  in¬ 
dividual  fracture  mode  are  determined  also.  For  given  adhe rends ,  the 
reduction  of  E^  increases  the  coefficients  of  K^.  and  in  Eq.  9, 

which  influence  the  G  value  significantly.  The  G^.  is  found  to  increase 
rapidly  with  E]/E2’  GII’  °n  t*ie  ot^er  hand,  decreases  as  a  less  stiff 
adhesive  is  used.  The  rapid  increase  of  G^  and  G  corresponds  to  a  decrease 
of  and  in  the  interface  crack  problem — a  unique  phenomenon  not 

observed  in  the  homogeneous  crack  problem  in  general.  In  the  case  of  very 
large  E^/E^  »  adhesive  may  become  incompressible  with  Poisson's  ratio 

approaching  0.5.  The  incompressibility  of  the  soft  adhesive  is  not  considered 
in  the  present  study. 

4.3  Effects  of  Dissimilar  Adherends 

Effects  of  dissimilar  adherends  on  the  interface  crack  behavior  are 
studied  for  lap-shear  joints  with  the  same  geometry  and  crack  length  used 
in  the  previous  section.  The  adhesive  layer  and  the  lower  adherend 
have  the  following  elastic  properties:  E?  =  0.5  x  106  psi,  =  0.35  and 
E^  =  10  x  106  psi,  =  0.33.  The  problem  is  solved  by  considering  the 
joints  with  various  E^'s  ranging  from  E-|/E2  =  1  to  1000.  The  results 
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are  given  in  Figs.  7  and  8,  in  which  ,  K  ,  G^,  G^  and  G  are  related 
to  E^/Ei  in  semi-logarithmic  plots.  Adherends  with  dissimilar  properties 
are  found  to  have  significant  effects  on  the  crack-tip  response.  The 
opening-mode  stress  intensity  factor  increases  appreciably  but  the  shearing¬ 
mode  stress  intensity  remains  relatively  constant  as  the  stiffness  of  the 
upper  adherend  approaches  that  of  the  lower  one,  i.e.,  E^/E ^  ->-20.  As 
E^/E^  exceeds  20,  an  opposite  trend  of  changes  of  and  are 

observed;  K^.  remains  relatively  unchanged  and  increases  gradually. 

The  ratio  of  K^/K^  in  the  joint  with  dissimilar  adherends  is  always 
smaller  than,  or  at  most  equal  to,  that  of  a  joint  with  the  same  upper  and 
lower  adherends.  Effects  of  dissimilar  adherends  on  energy  release  rates 
are  shown  in  Fig.  8.  The  is  approximately  one  order  of  magnitude 

larger  than  G^  in  the  entire  range  of  ^^/E^  >  20.  The  total  energy 
release  rate  increases  rapidly  initially,  and  then  remains  relatively  constant 
with  the  change  of  E^/E^.  As  the  stiffness  of  the  upper  adherend  exceeds 
that  of  the  lower  one,  crack  resistance  of  the  joint  becomes  less  sensitive 
to  the  dissimilarity  of  the  two  adherends.  But  a  higher  fracture  resistance 
is  obtainable  as  the  upper  adherend  becomes  less  stiff. 

4.4  Effects  of  Adhesive-Layer  Thickness 

Effects  of  the  adhesive- layer  thickness  on  fracture  of  the  joint  have 
long  been  recognized.  Bascom,  et  al.,  for  example,  [13,14]  showed  experi¬ 
mentally  that  fracture  energy  release  rates  are  related  closely  to  the 
bond  thickness.  Extensive  analytical  studies  on  debonded  composite  lami¬ 
nates  were  conducted  by  Erdogan  [5,6,29].  Several  unique  features  of  the 
geometric  variable  on  the  composite  failure  behavior  were  revealed.  Wang, 
et  al.  [27]  studied  center-of-bond  cohesive  cracks  and  reported  fundamental 
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characteristics  of  thickness  effects  on  the  adhesive  joint  fracture.  In 
this  portion  of  the  study,  lap-shear  joints  with  a  geometric  configuration 
and  crack  length  identical  to  the  one  used  in  previous  sections  are  considered. 
The  joints  are  made  of  identical  aluminum  adherends  and  epoxy  adhesive  with 
elastic  constants:  E^  =  E^  =  10  x  106  psi,  =  0.33,  and  E^  =  0.5  x 

106  psi,  =  0.35,  respectively.  The  adhesive  thicknesses,  ranging  from 
t^/40  to  t^/5  (i.e.,  from  0.00125  in.  to  0.01  in.),  are  considered. 

Relationships  among  K  ,  K  j,  G^,  G^,  G  and  the  normalized  bond  thick¬ 
ness,  t2/t1  (with  t1  =  t^  =  0.05  in.)  are  shown  in  Figs.  9  and  10.  The 
is  found  to  be  insensitive  to  the  change  of  the  adhesive-layer  thickness, 
but  increases  rapidly  with  t^-  Quantitatively,  is  always  larger  than 
K  in  all  cases  studied.  This  is  particularly  true  for  the  joint  with  a 
thinner  adhesive,  in  which  has  a  value  about  one  order  of  magnitude  higher 
than  that  of  in  general.  The  Kj  clearly  dominates  the  fracture  of  the 
joint  with  a  thinner  bond,  but  the  shearing  mode  failure  becomes  increasing¬ 
ly  important  as  the  adhesive  thickness  increases.  Associated  energy  release 
rates  are  given  in  Fig.  10  for  joints  with  various  t2/t^'s.  The  decrease 
of  the  adhesive-layer  thickness  leads  to  an  increase  of  fracture  resistance 
in  the  joints  (i.e.,  a  reduction  in  the  total  energy  release  rate  and  the 
maximum  cleavage  stress) .  The  change  of  fracture  resistance  with  the  bond¬ 
line  thickness  in  a  lap-shear  joint  is  consistent  with  the  results  reported 
in  Ref.  [29].  However,  variations  of  interfacial  crack-tip  stress  intensity 
factors  and  associated  energy  release  rates  are  different  from  the  solutions 
obtained  for  center-of-bond  cohesive  crack  problems,  in  which  K’s  and  G's 
are  almost  independent  of  the  adhesive-layer  thickness.  This  situation  is 
expected  as  stress  singularities  are  different  in  the  two  cases. 
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4.5  Effects  of  Crack  Length 

Effects  of  the  interface  crack  length  are  of  significant  interest 
because  of  the  particular  joint  configuration  and  the  loading  condition 
applied.  The  geometry  of  the  lap-shear  joint  studied  in  this  section  is 
identical  to  that  in  the  previous  cases .  The  two  adherends  are  the  same 
aluminum,  and  the  adhesive  is  an  epoxy  resin  with  material  constants 
reported  previously.  Solutions  obtained  for  various  crack  lengths  are 
reported  in  Figs.  11  and  12,  in  which  K  and  G  are  related  to  a/L.  Both 
K  and  K.^  shown  in  Fig.  11  increase  rapidly  for  shorter  cracks  (e.g. , 
a/L  <  0.02),  and  change  almost  linearly  as  tt/L  becomes  longer.  The 
opening-mode  stress  intensity  factor  is  approximately  four  times  higher 
than  that  of  the  shearing-mode  in  the  entire  range  of  a/L  studied.  Effects 
of  the  crack  length  on  energy  release  rates  are  shown  in  Fig.  12.  Well 
defined  relationships  are  observable.  In  the  case  of  a  very  short  crack, 
i.e.,  a  -*  0  ,  there  exists  a  strong  interaction  between  the  crack  tip  and 
the  edges  of  the  joint.  Basic  nature  of  the  short  interfacial  crack  and 
its  interaction  with  geometric  boundaries  are  currently  under  investigation. 


16 


5.  SUMMARY  AND  CONCLUSIONS 


An  investigation  of  the  interfacial  crack  behavior  in  adhesively 
bonded  lap-shear  joints  is  presented.  The  method  of  analysis  is  formu¬ 
lated  on  the  basis  of  conservation  laws  in  elasticity  of  nonhomogeneous 
solids  and  fundamental  relationships  in  fracture  mechanics  of  interface 
cracks.  The  current  approach  provides  a  convenient  and  accurate  means  to 
examine  the  basic  nature  of  interface  cracks  in  adhesive  joints.  Fracture 
parameters  such  as  stress  intensity  factors  and  associated  energy  release 
rates  describing  crack-tip  deformation  and  stresses  are  determined.  Complex 
failure  modes  in  the  joints  are  studied.  Solutions  are  obtained  for  problems 
with  various  kinds  of  adhesive  and  adherends,  joint  configurations  and 
interfacial  crack  lengths.  Based  on  the  results  obtained,  the  following 
conclusions  may  be  drawn: 

1.  Adhesive  properties  have  significant  effects  on  the  interfacial 

crack  tip  response.  In  the  lap-shear  joint  with  a  given  adherend,  the 
opening-mode  stress  intensity  factor  is  always  larger  than  that  of  the 
shea ring- mode .  The  difference  between  K  and  increases  with 

gl 

Ef/Ei"  While  both  and  decrease  with  decreasing  ,  the  total 

energy  release  rate  and  that  associated  with  the  opening-mode  fracture 
increase  rapidly  -  a  phenomenon  unique  to  the  interface  crack  and  not 
observed  in  homogeneous  crack  problems  in  general. 

2.  The  change  of  two  identical  adherends  to  dissimilar  ones  influences 


appreciably  the  failure  behavior  of  the  joints.  K  increases  rapidly  but 
K  remains  relatively  constant  as  the  stiffness  of  the  upper  adherend 
approaches  that  of  the  lower  one.  An  opposite  trend  in  the  changes  of 


K  and  is  observed  as  Z^/Z^  >20.  As  Z^/Z^  exceeds  20,  crack 


resistance  of  the  joints  becomes  relatively  insensitive  to  the  dissimilarity 
of  the  two  adherends •  The  joint  has  a  higher  fracture  resistance  as  the 
upper  adherend  becomes  less  stiff  than  the  lower  one. 

3.  Effects  of  the  adhesive  layer  thickness  on  interfacial  fracture  are 

different  from  those  on  a  center-of-the-bond  cohesive  crack  in  lap-shear 
joints.  Kj  is  shown  to  be  insensitive  to  the  change  of  the  adhesive  thick¬ 
ness,  but  increases  appreciably.  The  decrease  of  the  adhesive  thick¬ 

ness  leads  to  an  increase  of  fracture  resistance  in  the  joint  due  to  reduc¬ 
tion  of  crack  extension  driving  force. 

4.  Increasing  the  interfacial  crack  length  increases  K^.  and 
simultaneously.  However,  stress  intensity  solutions  increase  more  rapidly 
for  shorter  cracks  than  for  longer  ones.  For  the  cases  of  a/L  >  0.02, 
approximately  linear  relationships  between  K's  and  a/L  are  observed. 

In  a  semi- logarithmic  plot,  and  G  vary  with  a/L  almost  linearly 

too,  suggesting  that  well  defined  relations  among  them  may  exist. 
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